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We measured the neutron decay lifetime 
by counting in-beam neutron decay recoil 
protons trapped in a quasi-Penning trap. 
The absolute neutron beam fluence was 
measured by capture in a thin LiF foil 
detector with known efficiency. The com- 
bination of these measurements gives the 
neutron lifetime: T^ = (886.8 ± 1.2 ± 3.2) s, 
where the first (second) uncertainty is sta- 
tistical (systematic) in nature. This is the 
most precise neutron lifetime determina- 
tion to date using an in-beam method. 
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1. Introduction and Discussion 

As the simplest semi-leptonic decay system, the free 
neutron plays a crucial role in understanding the 
physics of the weak interaction and testing the validity 
of the Standard Model. The current experimental uncer- 
tainty in the neutron lifetime dominates the uncertainty 
in calculating the primordial helium abundance of the 



universe with Big-Bang nucleosynthesis [1]. Currently, 
there are six experiments [2-7] that contribute to the 
neutron lifetime world average of \ = (885.7 ± 0.8) s 
[8]. In the four most precise measurements, ultra cold 
neutrons were confined in a bottle; the decay lifetime is 
determined by counting the neutrons that remain after 
some elapsed time, with a correction for competing 
neutron loss mechanisms. The other two (in-beam) 
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experiments [2, 6] measured the absolute specific activ- 
ity of a beam of cold neutrons by counting decay pro- 
tons. Given the very different systematic problems that 
the two classes of experiments encounter, a more pre- 
cise measurement of the lifetime using the in-beam 
technique not only reduces the overall uncertainty of \ 
but also provides a strong check on the robustness of 
the central value. 

We measured the neutron lifetime at the cold neutron 
beam NG6 at the National Institute of Standards and 
Technology (NIST) Center for Neutron Research, using 
the quasi-Penning trap method first proposed by Byrne 
et al. This method is described in detail in previous 
publications [9]. Figure 1 shows a sketch of the exper- 
imental configuration. A proton trap of length L inter- 
cepts the entire width of the neutron beam. Neutron 
decay is observed by trapping and counting decay pro- 
tons within the trap with an efliciency £p. The neutron 
beam is characterized by a velocity dependent fluence 
rate /(v). The rate A^ at which decay protons are detect- 
ed is proportional to the mean number of neutrons 
inside the trap volume 



A^„ 



v„£da/(v)- 



(2) 



A^. 



£„L 



= ^Jda/(v)-, 



(1) 



where A is the beam cross sectional area. After leaving 
the trap, the neutron beam passes through a thin foil of 
^LiF. The probability for absorbing a neutron in the foil 
through the *Li(«, /)''He reaction is inversely propor- 
tional to the neutron velocity v. The reaction products, 
alphas or tritons, are counted by a set of four silicon 
surface barrier detectors in a well-characterized geom- 
etry. We define the efficiency for the neutron detector, 
£„, as the ratio of the reaction product rate to the neu- 
tron rate incident on the deposit for neutrons with ther- 
mal velocity v^ = 2200 m/s. The corresponding effi- 
ciency for neutrons of other velocities is SoVjv. 
Therefore, the net reaction product count rate N^ is 



The integrals in Eq. (1) and Eq. (2) are identical; the 
velocity dependence of the neutron detector efficiency 
compensates for the fact that the faster neutrons in the 
beam spend less time in the decay volume. This cancel- 
lation is exact except for a correction due to the finite 
thickness of the *LiF foil (+5.4 s), and we obtain the 
neutron lifetime \ from the experimental quantities 
NJNp, £„, £p, and L. 

The proton trap was composed of sixteen annular 
electrodes, each 18.6 mm long with inner diameter of 
26.0 mm, cut from fused quartz and coated with a thin 
layer of gold. Adjacent segments were separated by 3 
mm-thick insulating spacers of uncoated fused quartz. 
The dimensions of each electrode and spacer were 
measured to a precision of ±5 |am using a coordinate 
measuring machine at NIST. Changes in the dimension 
due to thermal contraction are below the 10"* level for 
fused quartz. The trap resided in a 4.6 T magnetic field, 
and the vacuum in the trap was maintained below 10"' 
mbar 

In trapping mode, the three upstream electrodes (the 
"door") were held at +800 V, and a variable number of 
adjacent electrodes (the "trap") were held at ground 
potential. The downstream three adjacent electrodes 
(the "mirror") were held at +800 V. We varied the trap 
length from 3 to 10 grounded electrodes. When a neu- 
tron decayed inside the trap, the decay proton was 
trapped radially by the magnetic field and axially by the 
electrostatic potential in the door and mirror After 
some trapping period, typically 10 ms, the trapped pro- 
tons were counted. In counting mode, the door elec- 
trodes were lowered to ground potential, and a small 
ramped potential was applied to the trap electrodes to 
assist slower protons out the door The protons were 
then guided by a 9.5° bend in the magnetic field to the 
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Fig. 1. Schematic representation of the experimental method (not to scale). 
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proton detector held at a high negative potential (-27.5 
kV to -32.5 kV). After the door was open for 76 |as, a 
time sufiicient to allow all protons to exit the trap, the 
mirror was also lowered to ground potential. This pre- 
vented negatively charged particles, which may con- 
tribute to trap instability, from accumulating in any por- 
tion of the trap. That state was maintained for 33 |as, 
after which the door and mirror electrodes were raised 
again to +800 V, reverting to the trapping mode, and 
another trapping cycle began. Since the detector need- 
ed to be enabled only during extraction, the counting 
background was reduced by the ratio of the trapping 
time to the extraction time (typically a factor of 125). 
Figure 2 shows a plot of proton detection time for a typ- 
ical run. 

Protons bom in the grounded electrode region inside 
the trap were trapped with 100 % efficiency. However 
protons that were bom near the door and mirror (the 
"end regions"), where the electrostatic potential is ele- 
vated, were not all trapped. A proton bom in the end 
region was trapped if its initial (at birth) sum of electro- 
static potential energy and axial kinetic energy was less 
than the maximum end potential. This complication 
caused the effective length L of the trap to be difficult 
to determine precisely. To avoid this complication we 
varied the trap length. The shape of the electrostatic 



potential near the door and mirror was the same for all 
traps with 3-10 grounded electrodes, so the effective 
length of the end regions, while unknown, was to good 
approximation constant. The length of the trap can then 
be written L = nl + L^^^, where n is the number of 
grounded electrodes and / is the physical length of one 
electrode plus an adjacent spacer. Zj„j is an effective 
length of the two end regions; it is proportional to the 
physical length of the end regions and the probability 
that protons bom there will be trapped. From Eqs. (1) 
and (2) we see that the ratio of proton counting rate to 
alpha counting rate is then 






p 

£nV„ 



^ 



("l+knil 



(3) 



We fit Np/Na ^s a fiinction of « to a straight line and 
determine T„ from the slope, so there is no need to know 
the value of L^^, provided that it was the same for all 
trap lengths. Figure 3 shows raw data from a typical 
run, proton count rate vs trap length n. 

Because of the symmetry in the Penning trap's 
design, L^^ was approximately equal for all trap lengths 
that we used. However there were three trap-length- 
dependent effects that broke the symmetry: the gradient 
in the axial magnetic field (the dominant effect), the 
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Fig. 2. Atypical plot of proton detection time after gating on the detector (t= 0). Regions 1 and 111 are used to sub- 
tract background from the peak region II. After correcting for deadtime in the time-to-digital converter, the result- 
ing peak area gives the proton rate N„. 
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Fig. 3. Typical raw proton count rate N^ vs trap length data, fit to a 
straight line (top), and residuals (bottom). These data have not yet 
been corrected for nonlinearities. 



divergence of the neutron beam (a minor effect), and 
the sHght variation in geometry of the individual elec- 
trodes (a very small effect). Each of these caused L^^ to 
vary slightly with trap length. A detailed Monte Carlo 
simulation of the experiment, based on the measured 
and calculated magnetic and electric field inside the 
trap, was developed in order to correct for these trap 
nonlinearities. It gave a trap-length dependent correc- 
tion that lowered the lifetime by 5.3 s. 

A variety of surface barrier (SB) and passivated ion- 
implanted planar silicon (PIPS) detectors were used to 
count the protons. The proton detectors were large 
enough so that all protons produced by neutron decay 
in the collimated beam, defined by CI and C2, would 
strike the 19.7-mm diameter active region after the trap 
was opened. The detector was optically aligned to the 
magnetic field axis, and the alignment was verified by 
scanning with a low energy electron source at the trap's 
center and with actual neutron decay protons. When a 
proton hit the active region of the detector, the efficien- 
cy for proton detection was less than 100 % due to pro- 
ton backscattering — a proton can backscatter without 
depositing enough energy in the active region to be 
counted above threshold. This is complicated by the 
fact that an uncounted, backscattered proton has some 
probability to be reflected by the electrostatic accelera- 
tion field back into the detector and counted. 

To determine the proton detection efficiency, we ran 
the experiment with a variety of detectors with differ- 
ent dead layer thicknesses and different acceleration 
potentials. The fraction of protons that backscatter were 
calculated using the SRIM 2003 Monte Carlo program 
[10]. We made an extrapolation of the measured neu- 
tron lifetime to zero backscatter fraction (see Fig. 4). 

The neutron detector target was a thin (0.34 mm), 
50-mm-diameter single crystal wafer of silicon coated 
with a 38 mm diameter deposit of ^LiF, fabricated at the 
Institute for Reference Materials and Measurements 
(IRMM) in Geel, Belgium. The manufacture of 
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Fig. 4. A linear fit of the measured neutron lifetime at varying val- 
ues of the detector backscattering fraction. The extrapolation to zero 
backscattering gives the free neutron lifetime. 

deposits and characterization of the ^LiF areal density 
were exhaustively detailed in measurements performed 
over several years [11]. The average areal density was 
p = (39.30 ± 0.10) |ag/cm^. The a particles and tritons 
produced by the ^Li(«, ;)'*He reactions were detected by 
four surface barrier detectors, each with a well-defined 
and carefully measured solid angle. The geometry was 
chosen to make the solid angle subtended by the detec- 
tors insensitive to first order in the source position. The 
parameter £q gives the ratio of detected alphas/tritons to 
incident thermal neutrons. It was calculated using 

e,=^\\Q.ix,y)pix,y)eix,y)dxdy, (4) 

46 •'•' 

where cTq is the cross section at thermal (vq = 2200 m/s) 
velocity, i2{x,y) is the detector solid angle, p{x,y) is 
the areal mass density of the deposit, and 6{x, y) is the 
areal distribution of the neutron intensity on the target. 
The ^Li thermal cross section is (941.0 ± 1.3) b [12]. It 
is important to note that we take the ENDF/B-6 1 <J 
uncertainty from the evaluation, not the expanded 
uncertainty, to be the most appropriate for use with this 
precision experiment. The neutron detector solid angle 
was measured in two independent ways: mechanical 
contact metrology and calibration with ^^'Pu alpha 
source of known absolute activity. These two methods 
agreed to within 0. 1 %. 

Proton and neutron counting data were collected for 
13 run series, each with a different proton detector and 
acceleration potential. The corrected value of the neu- 
tron lifetime for each series was calculated and plotted 
vs backscattering fraction, as shown in Fig. 4. A linear 
extrapolation to zero backscattering gave a result of 
T„= (886.8 ± 1.2[stat]±3.2[sys]) s. The summary of 
corrections and uncertainties is shown in Table 1. Our 
result will be improved by an independent absolute cal- 
ibration of the neutron counter, which would signifi- 
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Table 1. Summary of systematic corrections and uncertainties (in seconds) for the 
measured neutron lifetime 



Source 


Correction 


Uncertainty 


LiF foil areal density 




2.2 


Li cross section 




1.2 


Alpha detector solid angle 




1.0 


Sohd angle correction for beam distribution 


+1.5 


0.1 


LiF target thickness 


+5.4 


0.8 


^LiF distribution in target 


-1.7 


0.1 


Neutron losses in Si wafer 


+1.3 


0.5 


Neutron beam halo 


-1.0 


1.0 


Trap nonlinearity (Monte Carlo) 


-5.3 


0.8 


Proton backscatter calc. 




0.4 


Proton counting statistics 




1.2 


Neutron counting statistics 




0.1 



Total 



+0.2 



3.4 



cantly reduce the two largest systematic uncertainties, 
in the *LiF foil density and *Li cross section. A cryo- 
genic neutron radiometer capable of such a calibration 
at the 0.1 % level has recently been demonstrated [13], 
and we are pursuing this method further 
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